ii 











Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations l. Thesis and Dissertation Collection, all items 


1960 


The analysis of third-order servo with 
backlash and load inertia 


Sjamsoeddin, Abdoellah; Tanaka, Katsuji 


Monterey, California: U.S. Naval Postgraduate School 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
D U DLEY research materials and institutional publications created by the NPS community. 
qe Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
KNOX appointed — and published — scholarly author. 


LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


NPS ARCHIVE 


1960 
SJAMSOEDDIN, A. 





THE ANALYSIS OF THIRD-ORDER SERVO 
WITH BACKLASH AND LOAD INERTIA 
ABDOELLAH SJAMSOEDDIN 


and | 
KATSUJI TANAKA 


LIBRARY 
NAVAL POSTGRADUATE SCHOOL 


MONTEREY, CALIFORNIA 


us. 














UNITED STATES 
NAVAL POSTGRADUATE SCHOOL 





THESIS 


THE ANALYSIS OF THIRD-ORDER SERVO WITH 


BACKLASH AND LOAD INERTIA 


Abdoellah Sjamsoeddin 


and 


Katsuji Tanaka 





12ND P2238 (1-59) 





THE ANALYSIS OF THIRD-ORDER SERVO WITH 


BACKLASH AND LOAD INERTIA 


Жж 


Abdoellah Sjamsoeddin 
and 


Kat suf. Tanaka 





THE ANALYSIS OF THIRD-ORDER SERVO VITA 


BACKLASH AND LOAD INERTIA 


by 


Abdoellah Sjamsoeddin 
Captain, Indonesian Navy 
and 
Katsujs Tanaka 


Lieutenant (j.2.), Japanese Navy 


Submitted in partial fulfillment of the 


requirements for the degree of 


MASTER OF SCIENCE 
LN 


ELECTRICAL ENGINEERING 


United States Naval Postgraduate School 
Monterey, California 


1960 





THE ANALYSIS OF THIRD-ORDER SERVO WITH 
BACKLASH AND LOAD INERTIA 
by 
Abdoellah Sjamsoeddin 
and 


Katsuji Tanaka 


This work is accepted as fulfilling the 
thesis requirements for the degree of 
“MASTER OF SCIENCE 
IN 
ELECTRICAL ENGINEERING 
From the 


United States Naval Postgraduate School ድመ 








— 


7 


“ 


ҒАСЕ 


"OQ 
M 


Gacklach becomes 2n cs2íocering problem, fhe t.. i15 ordinarily arises 


Y 


іп the gear train connecting the output motoc ts» the isad, The ver- 


м 


formance kanges resulting from the addition of Sackizsh depend upon 
whethé. ihe major portion oft the output inertia is on the motor or the 
load side of the backlash. Wo matter where the secvomechanism error is 
messural, backlash contrikhutes either to the sreedy state error, or in- 
srabilit: bf the sy$tem, 
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l. Background of problem. 
In the study of backlash we have to consider two cases: 
a) Simple backlash theory, in which there is no inertia on 
load side. 
b) Complicated backlash i.e. inertia and friction on load 
Side. 
a. Simple backlash theory. 
A simple feedback control system is assumed as shown in Fig. l. 
Backlash occurs in gear train connecting the output motor to the load. 
The summer is assumed ideal here in its function of sending to the con- 
troller the difference between input signal and output of box containing 
gears with backlash. The ínput to the system is a shaft rotation, de- 


signated Ө in Fig. l. The output shaft is a rotation and designated 


R 


with 9c: The output of the gear box, whose output is 2O9,., is assumed 


С” 

ideal in all respect except for the presence of backlash. Because there 
is no inertia on load side the backlash is such a character as shown in 
Fig. 2. This figure shows that as the output shaft moves with finite 


velocity from a position ©, is zero to a position Ө. has some posi- 


C C 
tive value, say 9c1 greater than 1/2 the backlash angle, the output of 
the gear box 60% moves from zero to a certain positive value; then the 


E 
direction of rotation reverses so that e decreases through zero to 
some negative value €, greater than 1/2 the backlash angle. The value 
of ec goes from positive through zero to negative value and then goes 
back to zero again. ZN is the total backlash angle. There is a dead 
space and this means that the adjacent gear teeth must move through a 
finite distance before making contact with each other. As long as the 


operation is unidirection there is no trouble, but soon when it changes 


its direction, then backlash must be taken up each time. In manual 
1 





operation it is sometimes noisy and can cause instability in the system. 

b, Qualitative effect of inertia on load, 

Assumed we have a unity feedback system as shown in Fig. 1. 
The discussion of the physical nature of the system is as follows: 

In normal operation the motor may be driving in one direction and 
backlash is then taken up, so that all gears are in contact. If the load 
inertia is small and there is some load friction, the load shaft will stop 
almost immediately with the reversal of the motor rotation. The load shaft 
will remain stationary until the motor has moved through an angle suffi- 
cient to take up all backlash. On the other hand, if the load inertia is 
big enough, and the load friction is small, when the motor reverses the 
load will not reverse, but may continue to rotate in its original direc- 
tion with a constant velocity equal to its velocity at the instant the 
gear tooth contact was broken. If both the load inertia and friction are 

apprecisble, the load rotates in its original direction but with a con- 
stant deceleration. Under these conditions the motor shaft and load shaft 
both help to take up the backlash. 

с. Second order servo with backlash, 

The schematical diagram of a secondorder servo is shown in Fig. 
3. In the normal arrangement the load is geared to the motor through 
gear train с, and the load shaft is geared to a measuring device 
through gear train G5. There may be backlash in either of these gear 
trains or in both of them. For our purpose we consider the backlash in 
G, only. Consider the sketch of Fig. 4. Assuming a step displacement of 


magnitude OA and also assuming that the backlash happens to be all taken 


up in the proper direction, the trajectory starts at A. The trajectory 








Fig: 1 Block dia gram of а Simple Feedback 


control System. 


e 





Fig. 2. Relationship between angular position of two gears between 
which backlash exists. 


Load 





follows a linear path. At point D the velocities of the motor and 
load become different, mechanically contact is lost and separation 
occurs, providing the friction of load has a finite value and the coloumb 
friction of load can be neglected. The point of separation is determined 
by the load trajectory, when this trajectory has a common slope with the 
trajectory of the system. Combining these points together it forms a 
straight line and is called the "Separation dividing line." After 
point B the system operates in the backlash region. If the viscous 
friction of the load is negligible the trajectory of the load is a 
straight line with constant velocity. But when the friction is suffici- 
ent enough the load is drifting separately with a constant deceleration. 
The remainder of the system continues to drive in the reverse direction 
until the backlash has taken up and the system output completely cor- 
rected as to displacement. The separation dividing line determines the 
initial conditions for load velocity and displacement, which are the same 
value as the initial conditions for the motor velocity and displacement. 
The sketch of Fig. 5 shows the displacement of load as a function 
of time and also the motor displacement as a function of time. From 
these two curves the time at which the backlash is taken up can be deter- 
mined, Knowing this time the volocity and displacement of load at time 
of recombination can be determined. This is point C in the sketch of 
Fig. 4. At the same time point D is also located. This point indicates 
the velocity and displacement of the system without load. Each pair of 
corresponding points for load and system without load, velocity and dis- 
placement at recombination, in turn determine the velocity and displace- 
ment for the recombined system. This is obtained by applying the law of 
conservation of momentum. The point of recombination for the system is 


4 








Fig. 3. A second order servo with backlash 





e у” 


Trajectory of the combined system. Trajectory of the load only 





Fig. 4. 
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| 
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of the motor only. 
| 


Trajectory 


BOB’ = Separation dividing line, 


Phase plane trajectory of a second order servo with backlash. 





Load-displacement 


~ 


| 
! 
| amount of backlash 
| 


| time 


€ 


Fig. 5. A sketch of curves of displacement vs time. 





point E in Fig. 4. Trajectories for the recombined system start from 
E and the system retains its original characteristics until the next 
Separation dividing line is reached; separation and recombination will 
occur again as in the previous half cycle. If the final state of the 
system spiral into the origin the system is said to be stable. But when 
the final state of the system oscillates with a constant amplitude, the 
system is said to have a limit cycle. 
d. Higher order system. 
For higher order, because of the complex mathematics we have 
to break it up into several phase plane trajectories in order to analyze 
the system. For example: 
1, The plot of velocity vs. displacement (Fig. 6) 
2. The plot of acceleration vs. displacement (Fig. 7) 
3. The plot of torque vs. displacement (Fig. 8) 
The conditions at the point of separation is the same as the second 
order system, At the point of recombination, the law of conversation 
of mementum is still applicable, but we also have to consider the torque 
equation of the motor, in order to determine the acceleration. Knowing 
the acceleration ànd velocity the point of recombínation is located. The 
trajectory of the combíned system start again until the separation divid- 
ing line is reached. The picture repeats again as in the previous half 
cycle. 

Note here, that for the determination of the load displacement and 
velocity ít is not necessary to go back to the time function. In order 
to know the value of load displacement and motor displacement at which the 

backlash is taken up, just plot the load displacement versus motor displace- 


ment as shown in Fig. 19. Using this value we go back to Fig. 6, where 


8 





the velocity of motor is plotted against the motor displacement. For 
the calculation of velocity at recombination, the law of conservation 
of momentum is applicable as the case in the second order system, At 
the recombination point the torque equation of the motor is true, since 
the load has no torque at all and also providing the current of the 
motor is a continuous function of time during the separation. From 
Fig. 8, we can determine the value of the torque at the time of re- 
combination, Knowing the torque and the velocity, we can evaluate the 


value of acceleration at this point. 





velocity 


| 


و 
displacement‏ 


Fig. 6. A sketch of velocity versus displacement. 


10 





o—— == y 


displacement 


Figo 7. А sketch of acceleration vs, displacement of 
the combined system, 


11 





torque 


time 


Fig. 8. A sketch of torque vs. displacement of the motor. 





load displacement 


| 


time 


Fig. 9. A sketch of load displacement vs. motor 
displacement 
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2. Review of Second Order System with Backlash and Load Inertia. 
For our particular problem assume the gear ratio is equal to one, 
the equation for the combined system becomes: 
ev + Ji ) 8. + (5 + fj) 9 + KO, = КӨ (1) 


For the system with load removed: 


79 « £8. « K5. = KO (2) 


12, И ey yo (3) 
Given: Jy + J, = Ly + f, = 0.40 K = 1.0 
JM = 0.50 fy = 0.30 @, = 1.0 

J, = 0.50 f£, = 0.10 A = GB. OT 


Equation (1) becomes: 


16. + 0.406, + ec = 1 (4) 


To obtain isoclines for the combined system: 
6 1-9 
ec ቆሩ 
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Computed values of 9c and e for various values of N, are listed 


in Table one, 


The system as open loop with load separated: 
0.509, + 0. 300 = 1- ө, (6) 


At the instant of separation: 





At the instant of separation: 


E то 
ርዕ E мо ፡ oe 


Hence the equation (7) becomes: 


N, + 0.60 = ee 9 uo? (7a) 
9 
Mo 


where N = 


15 





Mo 


I) 


Mo N 0.60 


+ 


Upon separation of the load, equation (3) becomes: 


0.500, + 0.109, = 0 
N = 9, = i 0.20 
3 P. 
E 
Letting N, = N, = -0.20, equation (8) becomes: 
x a> SEO 

1 - uo 
hence Oo - (1- 0.20) 5, 
arc. tan. Mo = 78.7? 

l - Buo 

f - 101.3° 


Solve for equation )9) as a function of time we obtain: 


2 - 
0.50 ( s 91 - sd, OL ) FF OTO se, 0 
5 ( ›58 + .1) 61፲ = ( .5s * .D8, + 
where Em = Өм em 
Equation (13) becomes: 
1 Ф 
9 —— S + 1 6 
L 5 Мо s (ag T Mo 


Lo 


(8) 


(9) 


(10) 


(11) 


0 


) (12) 


A3 


r (13) 


(14) 





substitute the value of equation (11) into equation (14) it becomes: 


Ф 
ዘ 
pa 
l 
кә 
CD 
+ 
Фф 


(15) 


= - 6720, (16) 


Oy 
(o 
> 
e 
о 
Ф 
1 
ም 
кеі 
| 
ኮን 
Со 
St 
=> 
о 
+ 
in 
O 
& 
nn 
кеі 


Values of d versus time (t) are tabulated in Table 1. 
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TABLE 1, 


2t 





t (Ее ss) 51 ECE) -2 91 

2 JOBS 196 212 . 004 

‚077 49285 - d . 185 

‚ б "LS » 965 -.2 4205 
ED , 194 ‚905 -.2 2405 
1.4 244 RS = 1.04 
158 ‚302 1:252 -.2 1532 
22 ‚ 356 1.78 2422 1.38 
2.6 ‚405 24205 -.2 14825 
ED ‚450 2,25 -.2 205 
Bn . 493 2.465 -.2 23265 
3.8 332 2.66 -.2 2,46 
4.0 . 920 2.315 መ 24:55 
Fig. 10 is a graph of equations (16) and (17) By appropriate ordinate 


Scaling. recombination times were computed as lísted in table number four; for 
corresponding valucs of 6 


Values of eu 8 


M’ 


I 


9, for dividing lines 


are listed in Table four. 


Substitute the value of equation (15) into equation (6) and solve for 94 
as a function of time: 


con 1.00 


(17) 


.6t 


Ө, = - 16 Out + 111.455 Du + 13,345 0,26 - 125 9,01 


М 


Values of Ou and e versus tíme (t) are tabulated in Table three. 
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TABLE 2 


ER nn mmm mm emm rm am ae a nm rm rr ERR mR Rr mmm сс ሙ፦ 


t ы 22 -125, 613,345 e ЛАН Əy -1 
> ‚9608 .887 -120 11.82 -3.2 111.455 075 
4 ,923 ‚787 -115.15 10.5 -6.4 111.455 „105 
.6 .887 ,698 -110.9 9.305 -9.6 111.455 .260 
8 .8525 „620 -106.5 8.28 -12.8 111.455 .430 
1.0 ,819 ,550 О 1-55 -16 111.455 .70 
12 ‚787 „487 - 98.25 0 -19.2 111.455  .505 
1.4 ,756 „432 = 94,5, (ӨШСЕ :22:4 . 1109089889 2305 
1.6 .726 ‚384 - 90.85 5.12 -25.6 111.455  .125 
1.8 ‚698 ‚ 340 - 87.25 4.545 -28.8 111.455 -.05 
2.0 ‚670 ,302 - 83.8 4.035 -32,0 111.495 -.310 
02 ‚644 ‚268 - 80.45 3.58 -35.2 1118455 615 
2.4 ,620 ,237 E 211 -38.4 111.455 -1.285 
2.6 ,595 ‚210 - 74.45 2,805 -41.6 111.455 -1.790 
2.8 ‚572 ‚187 - 22 029 72256 -44.8 111.455 -2, 355 
3.0 aen ,166 T 72022 -48.0 111.455 -3.125 
2 „528 „147 - 66.0 1.961 -51.2 111.455 -3.784 
3.4 ‚507 ‚130 2650 lke -54.4 111.455 -4.61 
306 „487 ‚116 - 60.95 1.55 -57.6 111.455 -5.545 
3.8 ‚468 ‚102 - 58.5 1.361 -60.8 111.455 -6.484 
4.0 „450 „091 - 56:299 9 192105 -64.0 111.455 -7,532 
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To satisfy the principle of conservation of momentum, applied the next 


equation: 


(CUL EJ.) e s E (18) 


From equation (18) 8,, апа Ө are known, hence we can solve this equa- 
tion. The values of ®, upon recombination corresponding to values of 
6 upon separation aré listed in Table 4 


Mo 5 ፇ7 
TABLE 4 

мо NC 

10 2274 
„8 ‚260 
‚6 ‚1216 
.4 029 
n 07125 
„1 -,0927 
. 08 -,0925 
. 06 -.1032 
. 04 -,0656 
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Fig. tl is the phase plane presentation of On pu 91» and 9, 
as a combined system and when operating in the backissh region. A step 
input of 0.80 is seen to separate the load from tle system at A, recombina- 
tion takes place when the load is at B, and the baiancing of momentum 
between load and the system without load, causes 9c; d of the re- 
combined system to originate from point C. Reseparation occurs at point 


D and the trajectory is seen to spiral into a limit cycle defined by 


points E, F, G and H. The magnitude of the limit cycle is 0.550. 
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3. Physical analysis of Third Order System with Backlash and Load 
Inertia. 


The backlash existing in the gear train between motor and load is 
placed inside the feedback loop as illustrated in Fig. 12. When the 


backlash is tsken up, the differential equations for the system are: 


UN (S, e J) 6, + Um (WE, + FDO, = To 19) 
١ kn (20) 
К.Е = 1, 8 ሩፍ | (21) 
E = КӨ - КӨ, (22) 


Substitute equation (22) into equation (21) it becomes: 


w 


K,K.®, - ጾ2ጾጋ*ር = De + PI, (23) 


Differentiate equation (19) we obtain: 


ም ፡- 


2 Ве 2 
a | | . 1 ገ ን ፦ y 
Te = 1/N(N Jy + 2296 + ]/ኮ (N ie + የ ነፀ (24) 


= 1 с 
апа LL, = 1/K, (T0) (25) 


2 or - 2 SS 2 4. 
| / 
1/K,N CH et 2226. + 1/K NCN E £, + P/K NN ل‎ 
2 к > 
P/K,R(N Ls s Oe - K K 8k - KK, 6, (26) 
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Rearranging equation (26) and multiplied both sides by к, it becomes; 


K, K/ K On 


2 вё 2 2 LIP 
1/N(N Ju + 3,30. + 1/N(N E + d + N РУ, + PJ de. 


2 e 
P/N(N fu + የ1 )9 + K, K,K,06 (27) 
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Fíg. 12. A block diagram of a third order servo with backlash 
and load inertia. 
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amplifier constant Volts/Volts 

motor torque constant ft-lbs/amp. 
error measurement constant Volts/rad. 
K /L = Volts/Volts-henry 


R/L = Armature resistance = Ohms/henry 
Armature inductance 


Armature resistance Ohms 


Armature inductance = henry 

Gear ration . = rad., motor/rad.output 
Inertia of system without load Sue ie ыы. 
Inertia of load ce" Tos በክ mat 


Friction of system without load 


Friction of load only 


Displacement of the combined system. 
Displacement of the load only. 


Input displacement. 
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Letting N = 1, equation (27) becomes: 
(28) 
This is a third order differential equation, hence we can not solve 
it by using phase plane technique. In order to be able to draw the phase 
plane trajectory we have to solve the above differential equation either 
by longhand or by using analog computer. Solving equation (28) as a fun- 


ction of time (t), 


= ገሪ ገረ ሪ 
Е- mE cu. PISE 4 CST CKKZ3) ^ Sen) 
$3. (3 4 POS (Pty )5 o 
E e Ors y CPU DIS кк] P 
ተ ር« ተ (ዓ/ፓ ተዉ ተ ር 
Su. (የ፲ ጋ СИЕЗ СТ” 


letting KK = 1 and De = 1/5, equation (29) becomes: 


Ө MIKK +S Itp) (РЕЈ. ц G0) 
P ST SF OVS € PS CPI ር ጋ 
(qu Cty 4 P3] 6.1 ЖО 
ei, ‹፥/፲+ጾ/5” + ( ፤/ፓጋደ A VI 


S[Cs+ a) CS+ b+ yc JCS +b -jed] (31) 


AD 


or 


[5+ С7+РЈ]6., + Qo 
(Ba) strebt ዓዓ. 


where: CAFR oz zb 
{РИ зарж 
ገሪ 


2 


acp A 


pn H 


1+ imitial con dition s = O , equation £3i) becomes 


b-jc2t 
a 1 5 | pt Е Bu J | 
Б yl acb* + € J) а? 2а) + б + с“ 2 )с (а - 5-999 
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EU - jc)t 


: | 
2jc(a - b + jc) | ШЕ) (32) 
б = 1/3 (33) 
(s + a) (s + b + jc) (s + b - je) 
Ө EL LU ES b - je)t " 5 - jc)t (е 
P Bee Ten: БҰТЫ "ТЕТІН 
a -2ab + b с 2jc (a - b - jc) 216( 8 -b+ jc) 


(3ሩ) 


From equations (32) and (34) _we can draw the phase plane trajectory 


of the velocity vs. displacement (see Fig: 13) 


en 


96 = (1/J)s (35) 


(s + a) (s +b + jc) (s + b- jc) 
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T a E оу СЕ e 3 
AM | el te 000 0080 a) 
ER, 2 ላ 
a - "3-5 -c) 2jc(a - b - jc) ¿jala - 5 - jc) GF 
(36) 
Fig. 14 is tne plot of acceleration vs., displacement of the combined 
systen, 
The equation for the drifting load is determined by che friction and 


inercia of tha load. 


9 + 2б, = 0 (37) 
. 

И И 7 (38) 
Ө Jt 


This slote determives the point of seperation. 


Solving equation (37) as a function of time (ti, it becomes; 


= e 

ө; аж 1. 61. + 1 u Lo (39) 

5 5(8 + а j 
where 

ú =; CA 

ን . በ. = Ө = initial condition of displacement 

Lo Co Mo d 
e & 2 = in nditior 'eloci 

Pu e: MN initial condition cf velocity 
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Fig. 15 is the plot of velocity vs. displacement of the motor and 
the plot cf displacemert of Load vs. displacement of che motor in the 
backlash region. 

From the curve of the load-displacement vs. the motor-displacemént 
we can determine the point of recombination for the motor. Knowing 
the displacement of the motor, the velocity of the motor at recombina- 
tion can also be determined. 


Applied the law of conservation of mementum: 
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From the above equation 9 and 9, are known quantities, so we 
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can evaluate the valuc of E 
For the determination of the initial condition of the velocity at 


the point of recombination we have to apply the torque-equation of the 


motor. 
Т. = JO + £9 (53) 


Multiply equation (51) by Jy and equation (50) by d and add 


it together we obtain: 
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Fig. 16 is the plot of the motor-torque vs. the motor-displacement. 
Knowing the motor-torque and the motor velocity at the point of re- 
combination, we can determine cthe motor-acceleration at that point from 
equation (53). 

From here on we wolve for equation (32) using the initial condition: 
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until it hit the point of separation again. 
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9c + 1.39. + 5. 4 9. = 1,0 (55; 
Where J = 1.0 f = 1,0 М = 1.0 
u = © б Гм = o 68 P = 1 o 0 
: 4 са 22 
21 ; ft 32 KK 1.0 
?n = ju 0 
ር 1 = 1 
G(s) = „ы o Mme —— AQ 
ЕП: + 1.88 + 1.0) s(s + .9 + 3.436) (s + .9 - .3436) 


Fig. 17 shews the root locus of tbe system to be investigated, To 
Solve equation (553, used an analog computer and the result is plotted 


in Fig: (18). 
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Ths load separated with the equation: 


dch ^ 328, >= 0 (56) 
The slope of the load is 
e —. 05.88 (57) 
5. 
L 
Equetioma (56) is plotted tm Fig. (18). 
The equation of the mctor after separation: 
Sech i ek E - 8 А ) 
ዓነ + 2.4879, * 1.1339, + 1.6670, = 1.667 (58) 


From the slope of the load we can find the initial condition for 


equation (58). 


Using these intrial 
the result is plotted ix 


From Fir. 


in our case these are: 


ይ 
“Lo 


« ቂህ." ከር 7 
“ 455 "ж о. 7 


conditions we solve for squatien (58) and 


Fig. (18). 


(23) for a backlash = .42 rad. 


94 = 2328 гай D, = -.07% rad/sec, 
е, ы 27688 rad е 6480 rad/sec. 
T, а 024500 ft-ibe. 





o 6(. 480) 


E 
450 - 


„258 =5 


O” 
L 


E ° 


+ 


From Fig, (28), for Ө 


C 

Initial conditions are: 
Eo ri 
Eo м 
Өс» = 


Using these as initial 


again. Valves of ፀሩ. 9 


lated om table (5), 


E ,4(,074) = 


o 192 


‚258  rad/sec, 


rade. 


= 258 „328 гаа. 


rad/sec, 9. = 
2328 rad. 
.258 rad/sec, 


2192 ТРК) 


conditions we solve for equation (55) 


and ቅ. at points of separation are tabu- 


with several values of backlash. 
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- „396 
„ 380 
-„390 


„376 


TABLE 5 


or 
330 


-, 342 


„510 
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Backlash 


1.0 


1.0 


1.0 





= መያ eee eee eee eee መ ë — ` emm. ` gn. ` gem, ` መጨ eee eee ee ele eee 


720 
1.240 
.768 
1,9200 
. 788 
1. 200 
. 800 
1.192 


„810 
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424 
67256 
(356 


-.316 


o 334 
20902 
. 280 
-.276 
„25% 
S0 
P230 
-.240 


220 


TABLE 5 continued 
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- А0 
‚226 
-4222 
2192 
-.204 
. 186 
-. 0 
„170 
~ NOS 


o 165 


Backlash 





TASLE 5 continued 
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ር ር 9c Backlash 


er Run e so ee ee A MA er Me kat ete > er nn Fi er er iA “> መሙ ል. eg ege а ጩሬ ee ш 


1,360 
„676 
1.260 
‚79% 
1.208 
„798 
1.164 
„825 
1.140 
2532 


1.120 
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„800 
1„175 
„846 
[155 
„852 
1132 


. 864 
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-,170 


„158 


-.324 
«264 
-. 244 
. 200 
-.194 
168 
-. 166 
‚144 


-, 142 


o 264 
-.234 
.172 
-.160 
.134 
-,140 
„120 
120 


9110 





felt کر‎ „ E AA A سے‎ 


1.360 
720 
222 
.820 


1.136 
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. 1200 
860 
1.120 
„910 
1.074 


„940 
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2424 
-.344 
„276 
-.220 
„170 
-.148 


. 124 


TABLE 5 continued 


. 298 
-. 180 
. 140 
-. 112 
4094 
-.064 
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Backlash 





Values of 61, 91, M? Өм» 9c 9. апа Т, at points of recombina- 


tion are tabulated on Table 6 for several values of backlash, 


TABLE 6 
ጠጦፒጠ 24" 350 —S ТЕ 5-2. ашы 
ES L 9, On Ty “ር "e Backlash 
1.804 „070 .804 -.842 -.615 -.477 -.233 1.0 
. 124 -.050 1. 124 . 904 082 2522 „212 1.0 
1.840 .076 6840 -.852 -.636 -.481 -.251 1.0 
1.798 2074 „898 -„800 -. 604 -.4504 -.244 „9 
„170 -.060 1.700 „790 „606 „%50 „246 መ 
1.794 „670 .894 -.788 -.590 -.445 -.234 .9 
1.790 . 080 990 -.750 -.596 -.418 -.262 48 
„86 -„072 „986 746 584 „123 „246 „8 
1.784 „082 1.084 - .692 - 1582 -„383 -„275 7 
. 206 - „056 „906 „664 960 ‚376 259 „7 
1.734 „072 1.034 - . 680 -.548 2:309 -.245 ei 
1.776 . 090 1.176 -.632 E -. 343 -.280 .6 
2258 -.072 5 828 „610 . 526 23992 ‚2562 „б 
1.670 068 1,070 -a316 550 “ 1964 - 92 „6 
„398 -„05& „998 2554 2440 . 3108 ‚1914 .6 
1.558 . 058 . 958 54299 -.422 -.2996 -.1824 .6 
1.762 . 104 1.262 -, 560 -.552 -.294 -.316 አጽን 
„274 - .080 2778 2544 . 494 22944 . 2584 > 
1.602 .082 1.102 -.512 -.448 -.2744 -.2286 25) 
„404 - 046 . 904 ‚488 .406 .2744 . 1864 3 
1.508 .066 1.008 - 484 - , 390 -.264 -.1788 > 
-472 - ,042 .972 . 460 . 370 2.2592 . 1626 19 
1.472 „055 972 -.476 - ,364 - „2636 -.1332 49 
„512 -,034 1.012 2448 . 334 «2552 . 1298 > 
1.430 . 048 „930 - ,434 - , 324 - „2412 -.131 229 
45236 -.032 1.036 5 o 328 o 2404 „1376 3 
1.414 . 041 „914 -.426 -0 310 -. 2364 -. 1208 3 


E = መመ == — መመ — ይመ == ma ጨ፦- me == amo omm Loss =:- есе = — зава шеге ::=ሙ cvm 21 u 8: ምጅ p =>” Ld edi (-.። E = ፎፎሙ መሙ EU um 9 


Figs. 18, 19, 20, 21, 22, 23, 24, 25, 26, and 27 show us the result of 
the combined trajectories with backlash from .1 till 1.Oted,In our case 
limit cycle occurs, Values of limit cycles versus backlash is tabulated 
on Table 7. From the table we come to the conclusion that, if we in- 
crease the amount of backlash the amplitude of the limit cycle will in- 
crease too, 


De 





TABLE 7 


For Gain-constant is 1,0 





Backlash Amplitude of limitcycle 
ol --- 
„2 ==. 
93 == 
ግ | „399 
(2 0435 
„б „9625 
of 2425 
„8 ‚760 
29 810 

1.0 » 840 
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4, Analog computer tecrriqu: 


For che combined system: 


o 





This {ís an equation of a summer with 4 inputs i.e. 


— 


9 6. ል:22 6... 


On С” C ር 


С u 
SÉ es 

C 
0, 2 - 

n 


Equations (60), (61) and 
Hence for the combined system 


addition 3 sign chargers, 


(62) are equations for an integrator. 


we need 1 summer, 3 integrators and in 


(60) 


(61) 


(62) 





Summary for the combined system, 

The output of amplifier 1 which in our case is ቅ. is fed for- 
ward into amplifier 2 ({ ап integrator). One output of this ampiifier 
is fed back into amplifier 1 through a sign changer amplifier 6, 
The other output is fed forward into another integrator i.e. amplifier 3 
Output #1 of amplifier 3 is fed back into amplifier 1, 
Output #2 is fed forward into integrator amplifier 4. The output of this 
is fed back through a sign changer into amplifier 1. The output of 
amplifier 2 and amplifier 4 are recorded into a X- Y plotter i.e. 
the output of amplifier 2 i.e. 6 is fed into the Y - axis and the out- 


C 


put of amplifier 4 i.e. ፀሩ is fed into the X - axis of the X-Y plot, 


So we have a recording of acceleration versus displacement of the com- 
bined system. 


ሠ 


The other recording is output of amplifier 7, which is Ө. 15 


C 
fcd into the Y - axis and the output of amplifier 4, which is 9. 
is fed into the X - axis of the X - Y plotter, From this recording 


we have the phase - plane trajectory of velocity versus displacement of 
the combined system 
The circuit diagram of the analog computer set-up of the combined 


system is shown in Fig. 30. 
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Pi 


g. 


30 
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A circuit diagram of computer 


set up 


en 


i 





1 for the combined system. 
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Calculation for the computer set up # 1l. 


Amplifier 1. 


Amplifier 2 


Amplifier 3. 


Amplifier Rm 


J 
W 


Gu: Bei 4 1,0 


ኒያ 
!! 


1.0 


ርፍ ጠር 
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1,0 


1,0 


1,0 x .3 


1,8 x 33 


1.0 


1.0 


1.0 


= . 300 


z «9240 





Amplifier 5. 


Amplifier 6, 


Amplifier 7. 


W = 


"10 


tt 


lí 


The components used ín this'computer set up 


К, = 1.0 M-Ohms 
R, = 1,0 Meg 

К, ፡ 93 Meg 

К, = ‚4 Meg 

Re = 1.0 Meg 

Re = 1.0 Meg 

Ry = 1.0 Meg 

Re = 1.0 Meg 

Ry = 1.0 Meg 
Rio = 1.0 Meg 


fl 


£2 


£3 


£4 


£5 


£6 


£7 
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20 


1.0 


Ze = 1.0 
55 = 1,0 
ео = 1.0 
# 1 are: 
= 1.0 Meg. 


= 1.0 microfarad 
= 1.0 microfarad 
= 1,0 microfarad 
= 1.0 Meg 
= 1.0 Meg 


= 1.0 Meg 





Тһе time scaling 1.0 sec, 


The magnitude scaling for the displacement: 50 Volts 1.0 rad. 
The magnitude scaling for the velocity 2 50 Volts = 1.0 rad/sec. 


1.0 -— 


The magnitude scaling for the acceleration: 50 Volts 


When the system separate, the equation of the load becomes: 








Ө, - - .86, (63) 

6 = _1 | A (64) 
қ 8 

ô = 1 ЕЯ (65) 


The equation of the motor after separation becomes: 








Oy = | «5470, - .587Ө - 1/50 Berg °м | (66) 
Өм = 1 || (67) 
8 
ሀሀ ህ 
6, = 1 [^ | (68) 
S ፦ 
9, = 1 id (69) 
S 
ge чар - ኤ 60, - .48 Л (70) 
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Equation (64) and (65) are the equations for integrators. 
For the load equation we need two integrators and one sign changer. 
For the motor we need 3 integrators, one summer and 2 sign changers. 
The summery of the computer set up # 2 is as follows: 
One of the output of amplifier 1, which ts е; ie fed back into amplifter 
l. The other one is fed forward into amplifier 2, This amplifier has an 
output of d (load displacement), The output of amplifier 3 is fed 
into the Y-axis of a X-Y plotter, while the X- axis gets its input 
from amplifier 7. So we have a recording of load displacement vs. 


motor displacement, 


Amplifier 4 has 4 inpute i.e. O, 0 б, and 9,. One of the 


R' L' 
output of amplifier 4 is fed back into amplifier 4. The other one ts 
fed forward into amplifier 5. There are two outputs of amplifier 5, 
one is fed back into amplifier 4 through a sign changer amplifier 8. 
The other one is fed forward into amplifier 6, Amplifier 9 has 2 inputs 
1,6. Ou and ó 
Recording points are; 
9, from amplifier 5 is fed into the Y - axis of a X - Y plotter. 
9, from amplifier 7 is fed into the X - axis of the X - Y plotter, 
Hence we have 4 phase plane trajectory of motor velocity versus motor dis- 
placement. 


Motor torque T, from amplifier 9 is fed into the Y - axis and 


M 
che motor displacement from amplifier 7 is fed into the X - 
axis of the X - Y plotter. 

From this we have the phase plane trajectory of motor torque versus motor 


displacement. The circuit diagram of the analog computer set up for the 


load and the motor after separation is shown in Fig. 31. 
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31 s círcuit diagram of computer 
after separation. 


set up # 2 for 


{Э 


| 
| - 
С ታተ 


the load and motor 
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Calculation for the computer set up # 2. 


Amplifier 1, 


Amplifier 2. 


Amplifier 3. 


Amplifier 4. 


Amplifier 5. 


Amplifier 6. 


Amplifier 7. 


Amplifier 8. 


10 


11 


8 


gi 


#5 8 


Ceg* Po 


1.0 


H 


1.0 


„800 


1.0 


2547 


2547 


„800 


1.800 


1.0 


ато" و‎ 


1.0 


„800 


.347 


.547 


.800 


1.8 x 5005 = .901 


db 
፪ 


AE 1.0 
ag = 1.0 
aio = 1.0 
811 = 1,0 





Amplifier 9, 


а „К 
Wa = 12 £9 = .480 aj) = „480 
12 
H e Bugs Fee = ,600 = ,600 
13 ” ——_ EE * 413 4 
R13 
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The component used in the computer set up # 2 are: 


Ri = 1.0 Meg Gei = 1.0 microfared 
R, = 1.0 Meg Ceo = 1.0 microfarad 
R4 = 1.0 Meg Res = 1.0 Meg 

R, = 1.0 Meg Ce, = 1.0 microfarad 
Re = 1.0 Meg 

Re = 1.0 Meg 

Ry = „5005 Meg 

Rg = 1.0 Meg Ces = 1.0 microfarad 
Ro = 1.0 Meg Lee z 1.0 microfarad 
Rig = 1.0 Meg Rez = 1.0 Meg 

R11 = 1,0 Meg Reg = 1.0 Meg 

Rio = 1.0 Meg Reg = 1.0 Meg 

R14 z 1.0 Meg 


The time scaling and magnitude scaling used in the computer set up 


# 2 are the same as used in the previous computer set up. 


15 





5. The result of investigation. 
We run experiments with different values of the system gain i.e. 
г) К = 2-95 


b) K 1.00 


H 


с) К а= 1,20 

The result is plotted in Fig. 32. This figure shows the graph of 
the amplitude of limit cycle versus the backlash with different values 
of K. The graph tells us, that when we increase the amount of backlash, 
the amplitude of the limit cycle will increase too, provided we keep the 
gain constant. Now if we increase the gain, letting the amount of back- 
lash constant, the amplitude of the limit cycle will increase. Notice 
that by decreasing the amount of the gain by .05 (i.e. from K = 1.0 
to К = .95), the curve moves more to the right as does it to the left 
if we increase the system gain from K= 1.0 to K s 1,20. The values 
of the limit cycle with different gain constant is tabulated on table 


8. 
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TABLE 8 


Backlash Amplitude of limit cycle in rad. 
in K=- „95 К = 1.00 K = 1.20 
rad. 
1.0 2598 „840 -- 
9 2540 „810 Ар 
08 0425 „760 -- 
of . 305 e „725 „970 
„б -- 5625 „900 
«ጋ == 435 . 790 
Set -- . 353 „671 
З == == (552 
„2 == == -- 
ol -- ን Wi 
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5. Discussion and Conclusion 

It is seen from the results of this study that backlash in a third 
order system may cause a limit cycle or may not cause a limit cycle 
depending on the amount of backlash and the system gain. With large 
amounts of backlash and low gain a limit cycle is unlikely. As the 
gain is increased it becomes more probable that a limit cycle will 
exist. 

If a limit cycle does exist the amplitude of the limit cycle depends 
on both the amount of backlash and the gain. At low gain the amplitude 
of the Limit cycle (in radíans) may be considerably less than the 
amount of backlash (in radians). Conversely, at high gain the amplitude 
of the limit cycle may exceed the amount of backlash. 

These conditions have been determined from study of a system with 
adjustable gain, but without variation of the inertia and friction, nor 
of the amounts of inertia and friction on each side of the backlash. 


These factors need further investigation. 
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